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ABSTRACT 


Seismic  research  projects  undertaken  by  three 
Research  Associates  supported  by  this  contract  are 
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INTRODUCTION 


During  the  first  year  of  this  program  3  research 
associates  (Dr.  E.  Husebye,  M.  Saito,  and  R.  Turpening) 
have  been  appointed.  They  have  been,  either  on  individual 
basis  or  in  cooporation  with  the  other  members  of  the 
Department  of  Geology  and  Geophysics ,  engaged  in  a 
variety  of  seismological  research  problems.  These  include 
both  theoretical  and  data  oriented  studies  and  utilize 
data  from  LASA  as  well  as  single  stations. 

The  Research  projects  that  have  been  completed 

during  the  year  or  still  in  progress  are  2  (1)  Synthesis 

of  dilatation  and  rotation  seismograms,  (2)  Particle 

\ 

motion-mode  filters,  (3)  Structure  of  the  earth's 
core  (4)  Partial  derivatives  of  phase  velocity  of 
surface  waves  with  respect  to  anisotropy  factors, 

(5)  Application  of  array  data  processing  techniques  to 
a  network  of  seismograph  stations,  and  (6)  Excitation  of 
free  oscillations  and  surface  waves  by  a  point  source 
in  a  vertically  heterogeneous  earth. 

These  studies  are  briefly  described  in  the  following 


pages  of  this  report. 
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REASEARCH  PROJECTS 

1«  Synthesis  of  Dilatation  and  Rotation  Seismograms 

(M.  Saito; 

Prom  an  array  of  identical  3-component  stations 
it  is  possible  to  compute  dilatation  and  rotation 
seismograms  by  carrying  out  the  spatial  differentiations 
numerically. 

Let  (x,  y,  z\  be  Cartesian  coordinates  with  z 
axis  positive  upward,  and  let  (u,v,w)  be  components  of 
displacement  in  (x,y,z)  direction.  The  purposes  of  this 
program  is  to  compute 


Qa. 

a  u 

BIT 

ST  — — - 

■t  — - 

dx. 

_  Tnr 

g 

dz 

ay 

from  ovserved  value  of  u  and  v.  In  these  equations  <5^, 
represents  area  dilatation  and  TS £  a  component  of  rotation 
around  the  z  axis.  6^  and  are  unaltered  by  a 

rotation  of  coordinates  around  a  vertical  axis. 
contains  only  SH  type  motions,  but  0a  has  P,  5V  and 
Rayleigh  wave  components.  Besides  0U  and  ^  have 
different  phase  relationships  to  vertical,  longitudinal 
and  transverse  displacements. 
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In  a  homogeneous  infinite  medium  a  plane  P  wave 
traveling  in  x-z  plane  can  be  written  as 

U  kx-  rz)J 

w -iv  -£z ] 

IT  --  0 


in  customary  notations.  Phase  angles  between  0^  and 
the  longitudinal  or  vertical  component  can  be  computed  from 

-0 


*• = iz 

KJ-  1  ~ 

Let  us  assume  V  >  o  so  that  seismic  waves  come 
from  below.  The  phase  anqle  between  0^  and  ZcT  is 
always  -  IT and  that  of  0^  and  20  is  depend 

ing  on  the  direction  of  propagation. 

For  SV  waves  in  a  homogeneous  infinite  medium 
we  have 


uz 


-  iis 


(SV  wave) 


Note  that  the  phase  angle  between  and  zJ"  differs 

by  if  from  that  of  P  waves. 
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SH  body  waves  and  the  Love  waves  have  the  same 
characteristics  concerning  ^  and  .  The  values 

of  the  phase  angle  are  given  in  Table  3. 

The  relationship  in  Tables  1  and  3  holds  even  when 
the  observation  is  made  on  a  free  surface.  Table  2,  on 
the  other  hand,  does  not  hold  true  for  an  incident  SV 
wave  beyond  the  critical  angle,  for  which  ill 

assume  either  ±  IT  or  °* 

Rayleigh  wave  displacements  at  a  free  surface  are 
given  by 

&  ~  6  fair-  Hr-  n/2 )  J 
v-  -  o 


Where  £  is  the  ellipticity  of  particle  orbits  at 
the  surface  and  depends  on  the  mode  of  propagation. 
Phase  relationships  for  this  case  are 


£  =  -H6r 


(Rayleigh  Wave) 


Note  that  the  phase  angle  iu  ^  7?  in 

contrast  with  that  of  body  waves,  and  that  it  depends  on 
whether  the  orbit  is  retrograde  (  £  '0  '  or  prograde 
(  6  y  0  )  as  well  as  the  direction  of  nropagation. 
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Since  observations  are  made  only  at  discrete  points 
in  space ,  differentiation  of  u  and  v  must  be  done  by 
interpolating  observed  values.  Let  u  and  v  je  inter¬ 
polated  by  a  set  of  functions  ft.  (?: .  !/)'• 

U  ^  HO,  V)  »  Z 

t)  -v  V-CXAP*  j  UJ) 

Then  dilatation  and  rotation  are  given  by 


A  simplest  estimate  of  coefficients  )  and 

'V~A. ir)  is  given  by  the  leact  square  method.  For  this 

purpose  it  is  convenient  to  adopt  matrix  formulation. 

Let  u*- **<$,*;&  and  be 
observed  values  of  u  and  v  at  j-th  station,  and  let 

U  and  V  be  column  matrices  made  from  u^  and  v^ ,  respect¬ 
ively.  Then  the  least  square  estimate  of  and  A~.  are 

given  by 
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where  A  and  B  are  column  matrices  made  of  Ci^  and 
The  matrix  is  defined  by 


1 


J 


Substituting  ..  and  B  we  find  that  estimated  value  of 
u  and  v  at  each  station  are  given  by 


and  that  *he  residuals  are  given  by 


i 


u 

v 
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Matrices  for  dilatation  and  rotation,  whose  components  are 
dilatations  and  rotations  estimated  at  stations,  are 
given  by 


where 

by 


and  are  derivatives  of 

z 

L 


defined 


*.  V  N  > 


Preliminary  computations  show  remarkable  difference 
between  dilatation  and  rotation  seismograms.  As  will  be 
expected  r.o  distortional  wave  arrives  at  P  arrival  time, 
and  both  distortional  and  dilatational  waves  arrive  at 
S  arrival  time.  Comparing  phase  angles  between  strain 
seismograms  and  vertical,  longitudinal  and  transverse 
component  seismograms,  one  can  identify  several  modes 
of  propagation  as  well  as  the  direction  of  propagation. 


Table  I.  Phase  angle  for  P  wave 


direction  of  propagation 


-X 

cA^6»/a) 

-*/z 

-viz 

OAflO  a/tf) 

-'"/z 

-'ir/it 

Table  2.  Phase  angle  for 

SV  wave 

direction  of  propagation 

+  x  -  x 

aA^(6a./iO  7 r/2 

#A<|  (^/uT)  TT/i 

Table  3.  Phase  angle  for  SH  (Love)  wave 

+  HC  ~ 

//')  -  ^  1  ^/z 

Table  4.  Phase  angle  for  Rayleigh  wave 

retrograde  {  6  ^0  )  prograde  (£>0) 

4  *  -  X  4  X  -  * 

177  ^2- 

om  (Ba/v)  )  Q  i<Tf  i'ir  0 
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2.  Particle  Motion  -  Mode  Filters  (R.  Turpening) 

Analog  -  optical  mode  filtering  schemes  have  been 
applied  to  seismic  data  to  enhance  rectilinear  motion  in 
the  presence  of  retrograde  elliptical  motion  (Turpening, 
1966).  The  system  involves  projecting  the  instantaneous 
particle  motion  ptierived  from  three  component  information) 


The  projection  P  (  tj  Oj  X  ^  that  occurs  along  a  certain 

finite  number  of  these  half  rays  is  systematically  (oC«4? A  ■f-  ft*  ) 

placed  on  film  in  a  variable  density  format.  Such  a  film 
is  then  processed  by  optical  diffraction.  (Jackson,  1965). 

Digital  equivalent  of  above  optical  filter  can  be 
constructed.  At  the  present,  the  algorithms  are  being 
prt  ..rammed  to  apply  these  filters  to  digital  data. 


I 


Filtta  JL 


Filt cr  2, 


The  analytical  expressions  representing  the  output 


of  these  two  are: 


ZLt\  t  Jm.  U  0Tfl(t)  +  rWt)]J 
-  ^/l  W[A (f) 4(n -  f>h(f/c«ti,yt<^ ] 

+  k  Ft  —  pI/Gg 

+i?FS(<4[pK(./)4*  jsyr<jf-pJju  yi“»vytJ/] 

—  filtered  projection  P(*,t) 

(3  ~  transform  variable  of  ** 

o C  *  ^^P7\  ■#'  Qx. 

n  *  k/Al 

*  *  IXAtf 

/s  •  ran s form  variable  of  t,  -  frequency 


where : 


-11- 


am- 

Kip,  Ky) 


Z(th  time  series  representing  vertical  ground  motion 
yj»  cosine  and  sine  transforms,  respectively  of  Z(f) 


time  series  representing  ground  motion  in  the 
direction  transverse  to  the  direction  of  noise 
propagation. 

cosine  and  sine  transforms,  respectively,  of  DTN(t) 

time  series  representing  ground  motion  in  the 
direction  of  noise  propagation 

cosine  and  sine  transforms,  respectively  of 
DN(t) 


(fSC^T  JU 


A 


f,  «= 

fa  m  -  «c 

ht* ' 

&C<g  •£**** 

#(*•)«  unit  step  function  at  a  ■  0. 

*  number  of  incremental  steps  of  fl 


i] 


These  expressiors  show  that  the  above  optical  filters 
operate  by  adding  certain  input  traces  that  have  been 
shifted  in  phase  by  ninty  degrees.  The  upper  signs 
are  taken  when  filter  1  is  considered  and  the  lower 
signs  for  filter  "  .  Thus  filter  1  will  pass 
rectilinear  motion,  and  prograde  motion.  Filter  2  will 
pass  rectilinear  and  retrograde  motion  while  attenua¬ 
ting  prograde  motion. 

Other  more  complicated  particle  motion  filters  are 
also  being  constructed  in  the  Fourier  plane.  If  the  filter 
exhibits  no  variation  in  the  f-  (3  direction  than  the 
filter  will  operate  strictly  on  a  “particle  motion  basis". 
However,  those  particle  motion  filters  thatepperate  by 
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frequency  filtering  using  judgements  on  particle  motion 
may  be  constructed  here  by  the  appropriate  variations  in 
both  the  f  and  Q  directions. 
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3.  Structures  of  the  Earth's  Core  <E.  Husebye)  ^ 

A  series  of  core  models  have  been  proposed  by 
Jeffreys,  Gutenberg,  Bolt,  Adams  and  Randall.  These 
models  agree  fairly  well  with  each  other,  except  for  a 
crucial  transition  zone  between  the  liquid  core  and  the 
inner  core.  Calculations  of  JT/j*  and  * Vi*  from  the 
velocity  distribution  of  the  above  models,  and  comparison 
of  these  values  with  observed  characteristics  of  core 
phases,  indicate  that  none  of  the  cere  models  are  quite 
satisfactory.  For  example,  the  precursors  implied  in 
Bolt's  and  Adams'  and  Randall's  models  have  much  smaller 
amplitudes  relative  to  the  PKIKP  phase  than  should  be 
expected. 

Preliminary  results  of  our  studies  includes  observa¬ 
tion  of  PKIKP  down  to  about  some  irregularities 

in  dT/dA  in  the  distance  interval  110"  -  140"  with  at 
least  one  caustic  near  A  **  121°  .  Precursors  have  been 

observed  between  A*  105"  -  142",  being  most  prominent 
between  105"  -  120"  and  130"  -  142".  These  precursors 
have  usually  no  sharp  onsets  nor  large  amplitudes  contrary 
to  the  normal  behavior  of  core  phases.  An  alternative 
explanation  of  these  phases  might  be  P  wave  leakage 
into  the  core  from  the  diffracted  P  wave  propagating  in 
the  mantle-core  boundary. 

In  the  crucial  140"  -  150"  distance  interval,  our 


n 
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observations  favor  a  continuous  wave  velocity  increase 
in  the  transition  zone  between  the  inner  and  outer  core, 
roughly  agreeing  with  that  of  Gutenberg.  The  amplitudes 
of  waves  traveling  in  the  outer  core  are  larger  than  those 
of  the  inner  core,  requiring  a  modification  of  the  velocity 

gradient  in  this  part  of  the  core. 

The  data  used  in  the  above  analysis  come  from 
North  American  WWSS  and  LRSM  stations.  At  present  we 
are  combining  this  network  of  stations  to  a  super 
array  to  compute  dT/dA  and  d2W  for  the  core  phase.. 
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4 .  Partial  Derivatives  of  the  Phase  Velocity  of 

Surface  Waves  with  Respect  to  Anisotropy  Factors. 

(M.  Saito) 

Perturbations  of  phase  velocities  due  to  small  changes 
in  earth  models  are  calculable  from  the  partial  derivative 
data  of  phase  velocity.  The  same  technique  can  be  applied 
to  compute  dispersion  curves  for  slightly  anisotropic 
models. 

A  transversely  isotropic  medium  which  is  symmetric 
around  x3  ■  z  axis  is  characteristized  by  five  elastic 
constants,  cll,  c33,  c44,  c12,  and  c13.  In  a  cylindrical 
coordinate  system  (x,S,  z)  with  z  axis  positive  upward 
the  stress-strain  relation  is  given  by 

i\r  ~  c//  * Yr  +  C/Z 

'  C'*  ehY-  C/Y  +  ‘Vj  ^ 

A*  '  %  e>r  +  %  V 
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In  a  homogeneous  infinite  medium  three  different 
types  of  plane  wave  solution  exist.  When  the  wave  is 
propagating  in  a  direction  perpendicular  to  the  ? 
axis,  three  velocities  are  givei.  by 


For  waves  along  the  z  axis,  they  are  given  by 


It  is  convenient  to  use  two  elastic  moduli  \  and 
and  anisotropy  factors  £  ,  (j>  and  defined  by 

A  +3/C  ~  Cy/ 

j'~  c>/ 

£  = 

cy/  -  Z 
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In  terms  of  the  eigenfunction  y,  as  defined  in  a 

VVo 

previous  paper  (Saito,  1967),  the  partial  derivative  of 
the  phase  velocity  of  Love  waves  is  given  by 


For  Rayleigh  waves  partial  derivatives  are  given  by 

Slight 


Jrjt 


(2) 


From  these  data  we  calculate  perturbation  of  C 
changes  in  £  and  ^  by 


due  to 


(3) 
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to  the  fi. st  order  approximation.  As  an  unperturbed  model 
we  may  choose  an  isotropic  model  if  the  perturbations  A f 
A*  and  A  3  are  small  enough. 

As  an  example  a  Gutenberg-Qullen  A  model  is  perturbed 
to  an  anisotropic  model,  "n  which  }  =  6=  1.1  and  *  =0.9 
between  depths  of  19  km  and  100  km.  The  results  are  given 
in  Table  1.  In  this  table  C  is  the  phase  velocity  for 
isotropic  model.  ACp  the  predicted  value  computed  from  (3), 
and  ac  is  the  exact  value  of  perturbation.  The  accuracy 
of  prediction  is  within  10%.  It  is  interesting  that  the 
predicted  value  A  cp  is  always  greater  than  the  exact  value. 

Table  T 


Love  Wave 

Rayleigh  Wave 

T (sec) 

AC 

ACp 

C 

AC 

AGp 

K*WjW 

20 

3.91858 

0.07224 

0.07778 

3.62059 

0.06671 

0.07382 

30 

4.12424 

0.11357 

0.11850 

3.86213 

0.05101 

0.06141 

40 

4.24760 

0.12252 

0.12956 

3.93965 

3.03180 

0.03895 

50 

4.31419 

0.11223 

0.12103 

3.96395 

0.02171 

0.02126 

60 

4.35778 

0.10288 

0.11093 

3.97672 

0.01555 

0.01832 

70 

4.39247 

0.09596 

0.10275 

3.98919 

0.01136 

0.01303 

80 

4.42356 

0.09071 

0.09629 

4.00466 

0.00833 

0.00936 

90 

4.45322 

0.0864S 

0.09108 

4.02403 

0.00618 

0.00678 

100 

4.48242 

0.08291 

0.08673 

4.04753 

0.00461 

0.00499 
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5,  Application  of  Array  Data  Processing  Techniques  to  a 

Network  of  Seismograph  Stations 

Seismic  Arrays  have  proved  to  be  very  effective  and 
useful  tools  in  seismology.  With  an  array  of  instruments 
and  appropriate  data  processing  techniques  significant 
gains  in  signal  to  noise  ratio  has  been  achieved.  Furthermore, 
with  Large  Apperature  Arrays  quantities  (such  as  dT/d&) 
can  be  measured  directly  and  interpreted  in  terms  of 
structure. 

Since  the  establishment  of  a  very  large  aperature 
array  is  expensive,  the  possibility  of  the  utilizing  a  net 
of  single  stations  as  a  super  large  aperture  seismic  array 
was  considered.  For  this  purpose  the  Fennoscandia  seismo¬ 
graph  network  was  chosen.  The  network  covers  a  large 
area  with  maximum  dimension  of  (800  km  and  crustal  structures 
vary  at  different  sites. 

For  the  evaluation  of  such  an  array  two  separate 
studies  were  conducted.  These  are  (1)  Similarity  of 
•wave  signals  frc  station  to  station  and  (2)  Signal-to- 
noise  ratio  improvement  by  simple  processing.  In  the  former, 
correlation  of  P-wave  were  computed  between  stations  of 
varying  distances.  It  is  found  that  signal  shapes  re¬ 
main  the  same  although  the  amplitudes  may  vary.  In  the 
second, the  array  was  phased  by  simple  delay-and-sum  operation 
and  significant  improvement  in  signal-to-noise  ratio 
were  observed  both  for  earthquakes  and  explosions. 
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6 .  Excitation  of  Free  Oscillations  and  Sui face  Waves 
by  a  Point  Source  in  a  Vertically 
Heterogeneous  Earth 

(M.  Saito) 

This  work  was  described  in  a  technical  report 
and  it  will  appear  as  an  article  in  the  J.  Geophys.  Res., 
%rithir  three  months.  The  following  is  the  abstract  of  the 
paper. 


Radiation  patterns  of  surface  waves  and  free 
oscillations  for  vertically  heterogeneous  elastic  media 
are  derived  for  arbitrary  sources  using  variational 
equations.  The  results  are  expressed  in  terms  of 
normal  mode  solutions  and  source  functions,  and  show 
that  additional  calculations  other  than  normal  mode  sol¬ 
utions  are  unnecessary  to  construct  radiation  patterns. 
Source  functions  for  a  single  force,  a  single  couple, 
and  double  couple  without  torque,  all  in  arbitrary 
directions,  are  derived. 
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